Formation and rapid evolution of domain structure at phase transitions in slightly 
inhomogeneous ferroelectrics and ferroelastics 
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We present the analytical study of stability loss and evolution of domain structure in inhomogeneous 
ferroelectric (ferroelastic) samples for exactly solvable models. The model assumes a short-circuited 
ferroelectric capacitor (free ferroelastic) with two regions with slightly different critical temperatures 
Td > T C 2, where T c i — T C 2 <C T c \, T C 2- We show that even a tiny inhomogeneity like 10 _5 K results in 
splitting the system into domains below the phase transition temperature. At T < T C 2 the domain 
width a is proportional to (T c i — T)/(T c i — T C 2) and quickly increases with lowering temperature. 
The minute inhomogeneities in T c may result from structural (growth) inhomogeneities, which are 
always present in real samples, and a similar role can be played by inevitable temperature gradients. 
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I. INTRODUCTION 

The idea that the phase transition in electroded short- 
circuited ferroelectric proceeds into homogeneous mon- 
odomain state E| is very well known. Similar result 
also applies to free ferroelastic crystals. However, it has 
never been observed. Surprisingly, both electroded fer- 
roelectrics and free ferroelastics do split into domains, 
although they should not. The present paper aims to 
answer why. 

It is generally assumed that in the finite non-electroded 
ferroelectric samples the domain structure appears in or- 
der to reduce the depolarizing electric field if there is a 
nonzero normal component of the polarization at the sur- 
face of the ferroelectrics |l],|2| (in complete analogy with 
ferromagnets ||), if the field cannot be reduced by ei- 
ther conduction (usually negligible in ferroelectrics at low 
temperatures) or charge accumulation from environment 
at the surface 0] . On the other hand, in inhomogeneous 
ferroelastics (e.g. films on a substrate, or inclusions of a 
new phase in a matrix) the elastic domain structure ac- 
companies the phase transition in order to minimize the 
strain energy, as is well understood in case of martensitic 
phase transformations || and epitaxial thin films 

In search for reasons of domain appearance in other- 
wise perfect electroded samples, which is not yet under- 
stood, we shall discuss (i) a second order ferroelectric 
phase transition in slightly inhomogeneous electroded 
sample and (ii) a second order ferroelastic phase transi- 
tion in slightly inhomogeneous free sample. These prob- 
lems have not been studied before. We consider an ex- 
actly solvable case of a system, which has two slightly 
different phase transition temperatures in its two parts. 
While the phase transition occurs in the "soft" part of 
the system, the "hard" part may effectively play a role 



of a "dead" layer [|l0| and trigger a formation of the 
domain structure in the "soft" part with fringe electric 
fields (stray stresses) penetrating the "hard" part. One 
has to check this possibility, but the behavior of the 
corresponding domain structure is expected to be un- 
usual: it should strongly depend on temperature since 
further cooling transforms the "hard" part into a "soft" 
one, while the first "soft" part becomes "harder" . Since 
the inhomogeneity is small, one might expect that the 
domains would quickly grow with lowering temperature. 
We indeed find a rapid growth of the domain width lin- 
early with temperature in the case of slightly inhomo- 
geneous short-circuited ferroelectric and free ferroelastic. 
This behavior is generic and does not depend on partic- 
ular model assumptions. Generally, the inhomogeneous 
ferroelectric systems pose various fundamental problems 
and currently attract a lot of attention. In particular, 
graded ferroelectric films and ferroelectric superlattices 
have been shown to have giant pyroelectric |1 1| and un- 
usual dielectric response [n2[. 



II. PHASE TRANSITIONS IN SLIGHTLY 
INHOMOGENEOUS FERROELECTRIC 

We shall first consider the case of slightly inhomoge- 
neous uniaxial ferroelectric in short-circuited capacitor 
that consists of two layers with slightly different critical 
temperatures, so that, for instance, a top part "softens" 
somewhat earlier than the bottom part does. We assume 
the easy axis z perpendicular to electrode plates, and 
make use of the Landau free energy functional for given 
potentials on electrodes ip a (zero in the present case) §| 
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where P z (Pj_ lis the polarization component along (per- 
pendicular to) the "soft" direction, index p = 1(2) marks 
the top (bottom) part of the film: 

A 1 =A, 0<z<k, 

A 2 = A + SA, -l 2 < z < 0. 

Here Al(2) = a (T ~ ^ci(2)) an d SA > (meaning 
T C 2 < T c i). The constant a = 1/Tq, where To ~ T at 
(T c ) for displacive (order-disorder) type ferroelectrics, 
T a t ~ 10 4 — 10 5 K is the characteristic atomic temper- 
ature. 

The equation of state is 5Flgd[P]/SP = E = — V<y9, 
where <p is the electrostatic potential, or in both parts of 
the film p= 1,2 : 



E z = -d z (p = A P P Z 
E ± = A ± P ± , 



bp; 



DV 2 ± P Z - gd 2 z P z , (2) 
(3) 



These equations should be solved together with the 
Maxwell equation, div(E + 4ttP) = 0, or 



(d 2 z + e a Vi) if = 4ird z P z , 



(4) 



where the dielectric constant in the plane of the film is 
e a = l + 4ir/A±. 

A. Loss of stability. 

We shall now find conditions for loss of stability of the 
paraclcctric phase close to T c \ with respect to mhomo- 
geneous polarization. The stability loss corresponds to 
appearance of a non-trivial solution to linearized equa- 
tions of equilibrium. Indeed, at the brink of instability 
the system is in neutral equilibrium, defined by linear 
terms. We are looking for a nontrivial solution in a form 
of the " polarization wave" , 
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We shall check later that the stability will be lost for 
the wave vector kl\ ^> 1 while the scale of change of P z 
with z is h so that V\P Z = k 2 P z > gd 2 P z ~ P z /lf, 
and the last term in the right-hand side of (^) should be 
dropped. Going over to Fourier harmonics indicated by 
the subscript k, we obtain for the Poisson equation: 



^ - e a k 2 Vk = 4wP' z 



(6) 



where the prime indicates derivative ( /' = df/dz, f" = 
d 2 f/dz 2 ). We can exclude P z k with the use of the lin- 
earized equation of state (0), which gives 



<p' k = (A p + Dk 2 )P zk . 



(7) 



Substituting 

this into @, we obtain (p'£ — y 'JJ — ' tp k = 0, where 
we have used \A + Dk 2 \ /An <C 1, which is always valid 
in ferroelectrics. We shall see momentarily that the non- 
trivial solution appears only when A\ + Dk 2 < 0, while 
A2 + Dk 2 > 0. The resulting system is 



Vik + x\k 2 Vik = 0, 



(8) 
(9) 



where 

xl 



-e a (A 1 + Dk 2 ) /4tt, xl = e« (A 2 + Dk 2 ) /4tt. The 
corresponding solutions satisfying the boundary condi- 
tions for electroded surfaces (<p = at z = li, fo) read 



<p lk = Fsinxik(z + h), 
ip 2 k = Gsinhx2fc(z - h). 



(10) 
(11) 



The boundary condition at the interface (z = 0) reads as 

f'lk _ ^2k 



A t + Dk 2 A 2 + Dk 2 ' 



(12) 



where we have used \Ax + Dk 2 1 /4n <C 1. We obtain from 
Eqs. (P)-(p^) the condition for a nontrivial solution 



Xit&nxikh =X2tanhx2fc^2, 



(13) 



which has a homogeneous solution k = and the in- 
homogeneous solution with k 7^ (E3), hence we have 
to determine which one is actually realized. The inho- 
mogeneous solution is easily found for X2M2 ^ 1, where 
tanh can be replaced by unity. Close to the transition 
X2/X1 ^ 1) an( i f ne solution is 
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when X2kl\ ^S> 1. This gives the condition of stability 
loss in the form \A\ = Dk 2 + - |gp . There is no solution 

for xl < 0- The minimal value of A for the nontrivial 
solution (the actual onset of instability, if the transition 
with k — does not occur earlier) is defined by 
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(16) 



where we have introduced the "atomic" size d a t ~ vf 
comparable to the lattice parameter. We obtain the cor- 
responding tiny shift in the critical temperature [see es- 
timates below Eq.@] T cl - T c ~ T d at / 'tj 2 h. Hence, 



2 



the system looses its stability with respect to an inho- 
mogeneous structure very quickly below the bulk tran- 
sition temperature. It is readily checked that the as- 
sumptions we used to obtain the solution are easily sat- 
isfied. Indeed, X2^h ^ 1 and x%kh 2> 1 both corre- 
spond to approximately the same condition when l\ ~ 
I2 ■ 5 A ^> — 4 yts ^r^ > meaning that the difference be- 

7r 1 / 2 e a '! 

tween T c should be larger than the shift of T c . 

Now we have to determine when the transition into 
inhomogeneous state occurs prior to a loss of stability 
with respect to a homogeneous polarization. The homo- 
geneous loss of stability corresponds to A = A h found 
from 



A h h + (A h + 6 A) h = 0. 



(17) 



For the inhomogeneous state to appear first, there must 
be A c > A h , or 8 A > 7r3/2 ('i+^) 4^1, This means that 

very tiny inhomogeneity in the sample is enough to split 
it into the domain structure, 



T c i — T C 2 — Tq 



(18) 



which, for a film 1 mm thick, is estimated as T a t fe ^5 

ta 1/2 (10 4 - 10 5 )1CT 7 K = (1(T 3 - 1(T 2 )K for displacive 
systems, and T c -^- < (l(T 5 - 1CT 4 )K for order- 
disorder systems. Certainly, such a small temperature 
and/or compositional inhomogeneity exists in all usual 
experiments. 



B. Domain structure at T c2 < T < T cl (A < 0, 

A + SA > 0). 

After stability loss the resulting "polarization wave" 
quickly develops into a domain structure, as we shall now 
demonstrate. The notion of the domain can be applied 
when the domain width a = ir/k c becomes comparable 
and larger than the domain wall thickness W ~ y/D/\A\. 
The relation W < a gives [see Eqs.([l5 16)] 
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1/2, 

e a h 



< 1. 



(19) 



This is the same tiny temperature interval where the 
present scenario unfolds, and system quickly goes over 
into domain state well above the lower transition tem- 
perature T C 2, if it is larger than the value defined by 
Eq- @. 

In the region below T c \ where the domain structure 
forms (as shown above, it occupies most of the tempera- 
ture interval T C \—T C 2), we can use the linearized equation 
of state 

E z = (A + 3BP^)(P Z - Poi) = -2A(P Z - Poi), (20) 
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FIG. 1. Schematic of the domain structure with the pe- 
riod 2a in inhomogeneous ferroelectric film of the thickness 
h + fa- Top and bottom layers have slightly different critical 
temperatures T cl > T c2 , T cl - T c2 < T clj T c2 . Slightly be- 
low T c i the top layer splits into domains with electric fringe 
field propagating into the bottom layer (fringe field shown as 
the hatched area in the top panel). The domains persist and 
evolve below T c2 when both layers exhibit a ferroelectric (or 
ferroelastic) transition (bottom panel). 



where |i"bi| = \J — A/ B is the spontaneous polarization 
in the top layer, which gives P zl = P m + 2^\E Z , P z2 = 

■j-E z , for the top and bottom layers, respectively. In 
this case the equation for the potential ip (0) reduces to 
a standard Laplace equation (e c d z + e a V^J tp — 0, with 
the boundary condition 



e c id z tpi - e C 2d z <p2 = AttPoi[x), 



(21) 



where e cl = 1 + 2tt/\A\, e c2 = 1 + A-k/A 2 . 

The spontaneous polarization in the top layer alter- 
nates from domain to domain as Pqi(x) — ±|Pqi| = 
±y/—A/B. We are looking for a solution in a form of 
a domain structure with a period T — 2a (Fig. [j]), 



P 01 (x) =J2 P oike lkx , tp(x) =Y^ k e lkx 



k 



(22) 



with k = k n — 2nn/T — irn/a, n — ±1, ±2, ... Going 
over to the Fourier harmonics, we can write the Laplace 
equations for both parts of the film as 



(23) 
(24) 



e c2V2fe - e a k 2 (fi2k = 0, 
with the boundary conditions at the interface z — 

<Plk = <P2k, eci^ifc - e c2P 2 k = 4vrF ife (25) 

The corresponding electrostatic (stray) field part of the 
energy is found as |l0| 



1 



dAu s ip (z = 0) , 



(26) 
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where a s is the density of bound charge at the interface 
, corresponding to only the spontaneous part of the po- 
larization Pqi(x), and integration goes over the area A 
between two parts of the film. We calculate this expres- 
sion by going over to Fourier expansion (|2^) and using 
the fact that in the present geometry a s (x) = — Pqi(x) 
(and, therefore, its Fourier component o~ s k = —Poik), 
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D, 



,1/2 



tj 2 coth 



tj 2 coth , 



(28) 



with k = nn/a, n = 1,2,..., similar to [[t3|. Note that 
here Poifc = 2|P i|/* 7m , n — 2j + 1, j = 0,1,... and 
zero otherwise. Adding the surface energy of the domain 
walls, we obtain the free energy of the domain pattern 



F 
A 



lih IQP^a 1 



(29) 



where D n — Dk n . Not very close to T c \ the argument 
of coth is . —kl\ > 1 even for the smallest k — n/a 

y £d 1 ~ ' 

what is checked by the subsequent result (Eq.^), so that 
. Minimizing the free energy, we 
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find the domain width 

" 2 1/2/ 1/2 



£ l/2 
( c2 



14C(3) 



AiZi 



1/2 



(30) 



where Ai = Ji/Pqi = da^A] 1 / 2 is the characteristic mi- 
croscopic length, and d at = ^-D 1 / 2 is comparable to 
a lattice spacing ( "atomic" length scale) . The expression 

@ is valid when x Jj^kl 1 > 1, or \A\ > 2d at / (nel /2 l^ , 
meaning that one has to be below T c by a tiny amount 
T c i — T > Tod at / (el^ 2 lij , estimated earlier. Note that 
close to T c i one obtains for the domain width 



a K 



n^J 2 
7V2C(3) 



d a th 



1/2 



(31) 



and this value does not depend on temperature. We shall 
formally refer to this result as the Kittel domain width. 
Incidentally, close to T c2 the domain width is a m 

1/2 



tM /2 4 / 2 

14C(3) 



-Aiii 

-1/4 



oc e^ 4 , which formally diverges oc 



(T — T e2 ) . However, in the vicinity of T C 2 the in- 
duced polarization in the formerly "hard" part has about 
the same value as the spontaneous polarization in the 



"soft" part, P Z 2 ~ Pqi- Then the equation of state 
in the bottom part becomes strongly non-linear, since 
the cubic term is much larger than the linear term, 
BP% w BP^ = AP m Pi AP z2 > A 2 P z2 , in the equation 
of state (since A ^> A 2 close to T c2 ), so the response of 
the bottom layer does not actually soften in this region. 
In this case our derivation does not apply, but it is prac- 
tically certain that the domain structure in the vicinity 
of T c2 would evolve continuously upon cooling, Fig. 0. 

C. Domain structure at low temperatures (T < T C 2, 
A< 0, A + 6A< 0). 

When the system is cooled to below the critical temper- 
ature T C 2, a spontaneous polarization IP02I = \J —A 2 /B 
also appears in the bottom layer. The domain structure 
simultaneously develops in the whole crystal with domain 
walls running parallel to the ferroelectric axis through the 
whole crystal (if they were discontinuous at the interface 
between the two parts of the crystal this would have cre- 
ated a large depolarizing electric field) . The electrostatic 
energy requires a solution of the same Laplace equations 
( p3| ) and (pif), only the boundary condition (|2q) would 
now read 



e cl¥>'lfc - £c2^ 2 fc 



47r(Foife - Po2fc), 



(32) 



where e cl ( 2 ) = 1 + 27r/| J 4 1 ( 2 )| ~ 27r/|A 1 ( 2 )|- Note that the 
density of the bound charge at the interface, correspond- 
ing to this discontinuity of spontaneous polarization, is 
now Ok = — (Poifc — Po2k)- Therefore, we immediately ob- 
tain for the total free energy of the structure, analogously 
to the previous case (p9f), 

F = P^h + P 2 2 A 2 l 2 
A a 



16(P i - P 02 ) 2 a 



y — -■ 

<T* (27 + 1)' 



% (2j + iy D 2j+1 



(33) 



where Aip) = d at ^J\Ay^\. Not very close to T c2 we 

would have . —kl 2 > 1 even for the smallest value of 

k = n/a which enables us to replace coth by unity. The 
minimum of the free energy F is achieved for the domain 
width 

1 



1 - P02/P01 
2 1/2 / 1/2 



+ e 



1/2 

<::2 



14C(3) 
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A2*2pP 



1/2 



(34) 



Close to the critical point T c2 the domain width formally 
behaves as a oc e]. 2 4 oc (T c2 — T) -1 / 4 , as found just above 
T c2 before. The same argument indicates though that 
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our derivation does not apply in this region, but non- 
linearity should not cause a substantial change in the 
domain structure. 

With lowering the temperature to the region where 
\A\ > SA, we will have P02/P01 = y/(A + 5 A) /A ps 1 + 
5A/2A, so that 1 - P02/-P01 ~ 2\A\/5A > 1 becomes 
a large prefactor. Note that in this region e c \ « e C 2 = 
2ir/\A\, Ai « A2 = dat and the domain width 
evolves as 



SA 



25/2^5/2^/2 

7((3) 



d a t(h + h) 



1/2 



(35) 



It becomes much larger than the Kittel width, 



2 3/2 



h+h 



1/2 



» 1, 



~Y — 1 rp - ^ -l, growing linearly with low- 
ering temperature (Fig. |^). For large periods of the 
domain structure Eq.(|35|) becomes inapplicable because 
the coth in the formula for Dk ( |28| ) cannot be replaced 
by unity, and this corresponds to \A\ ~ {5 A) 2 J 2 -. If we 
assume that the difference between the critical temper- 
atures in the both parts of the system is, for example, 
just T cl - T c2 = O.IK. Since \A\ = (T cl - T)/T and 
5 A = (T c i — T C 2)/Tq, we see that in 1mm thick film 
(I /d a t ~ 10 7 ) the expression for the domain structure 
period, Eq.([35|), is valid at least in the region 



T c i - T < 



h (T cl -T c2 ) 



(36) 



d a t Tq 

This interval is 1 — 10K for displacive and 50 — 100K for 
order-disorder systems. 

It follows from the qualitative analysis of the expres- 
sion for the electric energy, that the domain width a 
will keep growing with lowering temperature beyond this 
range to sizes much larger than the Kittel width (|3l|) , be- 
cause the system quickly moves into the region \A\ 3> 5 A. 
This result is rather natural, since in this limit the rel- 
ative difference between two parts of the system dimin- 
ishes, and the system approaches the limit of a uniform 
free sample, which transforms into a monodomain state 
(i.e. a = 00). 



III. INHOMOGENEOUS FERROELASTIC PHASE 
TRANSITIONS 

Very scenario applies in a case of slightly inhomo- 
geneous proper ferroelastic in spite of some differences 
in the mathematics. Consider a ferroelastic slab with 
slightly different phase transition temperatures, T c \ > 
T C 2, in its two parts of a comparable size. In such a 
situation, the "hard" part will play a role of a rigid sub- 
strate for the top "soft" part of the sample at tempera- 
tures slightly below T c \ , and the sample will split into do- 
mains. The emerging domain structure should strongly 



a/a 



cz cl 



FIG. 2. The domain width in slightly inhomogeneous fer- 
roelectric or ferroelastic in the units of clk, the Kittel width 
([jl|). a — ax when the domain structure sets in at T w T c i, 
and then it grows linearly with the temperature to large val- 
ues a S> ok . 

evolve with temperature, since the bottom part of the 
film would also become "soft" at T = T C 2 slightly below 

T C 1. 



A. Loss of stability (T w T cl ) 

We assume that the film is perpendicular to the z-axis, 
occupies the space — 1% < z < l±, and is characterized by 
the u xy (in-plane) component of the strain tensor as the 
order parameter. The "hard" shear modulus equals \i in 
both parts of the film, while the "soft" modulus corre- 
sponds to the u xy component of strain. We shall consider 
a situation when the system consists of two layers with 
slightly different critical temperatures, Fig. ||. Thus, the 
Landau thermodynamic potential has the form 



E 

P=l,2 



dV[2A p u 2 + 2D {Vu xy 



Bu ty 



(37) 



where A 1 = A = a(T-T cl ), A 2 = a(T-T c2 ) 
A + SA, with a, D, /i > positive constants, and SA > 
corresponding to T C 2 < T c \. Thus, the top layer of the 
system "softens" at T c \ (A — 0) while the other part of 
the system remains " hard" . Note that A designates now 
an elastic modulus and not reciprocal dielectric suscepti- 
bility as in the previously described case of ferroelectric. 
We designate the two parameters by the same letter to 
underline the similarities in the corresponding formulas. 
The equation of state in each part p of the system is 



\SF/Suik, i ^ k : 



2(Ap- 
2\xu xz , 



DV 2 )u xy 



2Bui 



(38) 
(39) 
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(b) 



FIG. 3. Schematic of the domain structure in inhomo- 
geneous ferroelastic film of the thickness h + I2 with "soft" 
in-plane strain u xy . Top and bottom layers have slightly dif- 
ferent critical temperatures T c \ > T C 2, T c i ~ T C 2 <C T c i,T C 2- 
(a) Slightly below T c i the top layer splits into domains with 
fringe elastic field near the interface z = (schematically 
shown by the dotted lines), (b) The domains persist and 
evolve below T C 2 when both layers exhibit a ferroelastic phase 
transition (bottom panel). 



where p = 1 (2) corresponds to the part < z < l\ 
(-l 2 <z<0). 

The treatment of the stability loss is analogous to that 
of Ref. |Q and we omit some preliminary discussion pre- 
sented in that earlier paper. To find the inhomogeneous 
part of the displacement vector u x (or, equivalently, u y ) 
at the phase transition one should satisfy the equations 
of local equilibrium, d<Ji k /dx k — 0, which in the present 
case read 



dcr xy _j_ dcrx 
dy 



dz 



0. 



(40) 



We shall use the Fourier expansion for the displacement 
vector 



u(y, z) — I u k (z) exp (iky) dk 



(41) 



and find the first appearance of the non-trivial solution 
for u for a given k wavevector. We then determine the 
k where the instability sets in first, and this will be the 
point of the stability loss of the symmetric phase. 

We obtain the following equations for the displacement 
with the use of Eqs.([4C 



d 2 u k A x 



k 1.2 



dz 2 
d 2 u k 
dz 2 



k*u k =0, 0<z<h; (42) 
-h < z < 0, (43) 
A+SA+Dk 2 . At the free 



^2k , 2 n 

k Uk = 



where A\k = A+Dk 2 and A^k 
surfaces (z — l\, —I2) the boundary condition reads a xz = 
0, which is equivalent to duk(z)/dz — 0. In addition, 



the displacement u k (z) and the stress a xz (z) should be 
continuous at the interface z — 0. 

Let us first consider the case of A lk < 0, A 2k > 0, 
which would correspond, as we will see shortly, to a loss 
of stability of the paraphase. The solution of Eqs. 
© is 

u k {z) = Fcosxik(h - z), < z < h; (44) 
u k (z) = G cosh X2k(z + l 2 ), ~h < z < 0, (45) 

where \\ = -Mk/p> = (-A-Dk 2 )/^ and \l = A 2k /^i = 
(A + SA + Dk 2 )/fi. 

The condition of existence of a non-trivial solution, 
which we obtain from the boundary conditions, looks 
exactly the same as in the above case of ferroelectrics, 
Eq.(|l3|).The subsequent analysis is also similar, and we 
obtain a condition for an existence of a nontrivial solution 



1^1 



2 

WI 2 



Dk 2 



(46) 



so the minimal value of \A\ = \A\ C , when the solution 
first appears, corresponds to 



\A\ C = 2Dk 



2 _ TTfldat 



(47) 
(48) 



where we have introduced the characteristic "atomic" 
1/2 

length scale d a t ~ (D/fi) ' , which is comparable to unit 
cell size. The corresponding shift of the critical temper- 
ature is very small. The coefficient a in (^) is fi/To, 
where To ~ T at in the case a displacive, and ~ T c in 
the case of order-disorder phase transition. Then, from 
Eq.©, 



T c — T c i 



KLld at 



T c i — Tq- 



Ttdn 



(49) 



which is practically the same estimate, as for ferro- 
electrics, with the same (by the order of magnitude) val- 
ues for the displacive and order-disorder phase transitions 
and the same condition for transition into inhomogeneous 
instead of a homogeneous state. 



B. Domain structure in the top layer at T C 2 < T < T c \ 

{A < 0, A + 5A > 0) 

We consider next the domain structure below T c \ in 
a state with the spontaneous strain u xy . One can apply 
the notion of the domain structure when the domain wall 
thickness is much smaller than the domain width. This 
condition is fulfilled just below T c \ by a tiny amount given 
by the same small parameter d a t/h as in the case of 



the ferroelectrics, Eq.(19). As in the previous case, one 
can apply the notion of the domain structure practically 



G 



in the whole interval T c % — T c2 , if the transition indeed 
proceeds into inhomogeneous state. Within this interval 
one can use the linearized equation of state for the top 
layer, obtained by expanding the free energy (^) about 
the spontaneous deformation, 



1 

xy 


xy 
2 



2M 1 (u xy -u° xy ), 

±(-A/B) 1/2 , 

-h < z < 0, 



< z < h 



u 



a- y = 2M 2 u xy , 



(50) 
(51) 
(52) 



where M\ = —2A > is -C n when the system is close 
to the transition ("soft" modulus), M 2 = A + 5A > 0. In 
both parts a xz is given by Eq.(f39|). 

We shall assume that all the domains have the same 
width |Q, which we will find by minimizing the sum of 
the elastic energy and the (surface) energy of the domain 
walls. We consider a stripe-like domain structure in the 
top layer with the spontaneous strain u xy (y,z) = ±uq 
with the period 2a. There would be no stresses in the 
free top layer if Uq = —A/B. We have to find the dis- 
placements u x (y, z) = u(y,z) appearing after the top 
layer experienced a phase transition. The equation of 
equilibrium ( |io|) takes the form 



„ r <9 2 u 8 2 u „ 
Mi— + M — = 2Afi— ^, 
oy z oz z oy 



, r „ "u a u 
M 2 —+n— = 0, 



'dy 2 



dz 



< z < h (53) 
h < z < 0. (54) 



Since the domain pattern is periodic, the elastic displace- 
ments may be represented as a Fourier series 



, (y, z) = y~] u k (z) exp (iky) , k = k n = — , (55) 



where n — ±1, ±2, . . . After solving the resulting system 
of ordinary differential equations with the above condi- 
tions one finds 



u k (z) = u R k 



r)2 coshrjik(z — li) 



771 smhrjikli coth r]2kl 2 + r\ 2 cosh 771 kl\ 



1 



at < z < h, and 



u k (z) 



uoRkVi cosh7/ 2 fc(z + l 2 ) 



r/i coshr] 2 kl 2 + 772 coth r/ikli sinh 772^^2 



(56) 



(57) 



at -l 2 < z < 0. where R k = 4:/ k 2 a, k = n (2r + 1) /a, 
r = 0, ±1,... and 77; = \J MjJJl, i = 1,2, The elastic 
energy is found with the use of the formula jl4| 

F -^~lf ^%dV = - J a xy u° xy dV (58) 

with the result for the elastic (stray) energy per unit area 
of the film: 



Fstray 



167717720 



E 



(2J 



(59) 



with Dj = T]iCothr] 2 kjl 2 + r\ 2 coth r\\kjl\. To find the 
equilibrium domain width, we have to add the energy of 
the domain walls 



Fdw/A = jih/a = /jTJoA^i/a, 



where 



7i 



Ai 



8V2D 1 / 2 \A\ 3/2 = 
8^D l / a \A\ 1/2 



3// 



= d at y/\A\/fl 



(60) 

(61) 
(62) 



and we have introduced the microscopic length scale Ai, 
with d a t = (D I ' [i) 1 ^ 2 is, once more, of the order of 
the unit cell size. Note that the actual domain walls 
exist only in the top layer, which underwent a ferroelastic 
transition, with the stray displacement field penetrating 
into the bottom "rigid" part of the sample. 

The equilibrium domain width is found from the total 
free energy 



F, 



stray 



Aiii 



(63) 



with the stray energy from (p9|). Assuming 
7T77i(2)^i(2)/ a ^ 1 (to be checked later), we replace all 
coth in (p9) by unity and easily obtain for the domain 
width 



'^Ai^i 771 + 772 



V 14C(3) mm 

Slightly below T c \ one has 771 <^ m 



(64) 



a = a>K 



I^Aih 1 



14C(3)»ji V 2 3 / 2 7C(3) 



dath, (65) 



the limiting value which does not depend on temperature 
close to transition. We shall formally call this a Kittel 
period for the elastic domain structure, and the system, 
as we have shown, looses stability and quickly sets in the 
domain structure with this period, which is independent 
of the temperature close to the phase transition. As in 
the case of ferroelectrics, the period of the domain struc- 
ture formally increases close to T C 2, but this conclusion 
is not reliable because the nonlinear effects in the former 
"hard" layer should be taken into account in this region. 



C. Ferroelastic domain structure at low 
temperatures (T < T c2 , A < 0, A + 5 A < 0) 

We consider next the domain structure not very close 
to the phase transition. There the domain wall width 
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is much smaller than the width of the domains and one 
can use the linearized equation of state in both top and 
bottom parts of the film, obtained by expanding the free 
energy (p7j) about the spontaneous deformation, 



a xy = 2M 1 {u xy - u° xy ), 
<j xy = 2M 2 (u xy -w° y ), 



< z < h 
-h < z < 0, 



(66) 
(67) 



where Mi = -2 A is < \i ("soft" modulus), M 2 = 
—2(A + 8 A) [note the change in the M 2 value below T c2 \. 
Similarly to previous case u xy = uq = ±(—A/B) 1 ^ 2 , 

w° xy ^w a = ±[-(A + SA)/B} 1/2 . 

We have to find the inhomogeneous displacements in 
the film u x (y, z) = u (y, zY For the film the equation of 
mechanical equilibrium ( [40|) takes the form 



d 2 u 



M 1 —+fi— = 2M 1 - 

oy z oz z ay 



Mo 



'dy< 



dz 2 



= 2M 2 



xy 



dy 



< z < h, 
-l 2 < z < 0. 



We look for a solution in the same periodic form 
earlier with the result 



Uk(z) = Rk 



(ito — wo) V2 cosh 771 k(z — h) 
rji smhrjikli cot}ir) 2 kl 2 + r\ 2 cosh 771 kl\ 



(68) 
(69) 
|) as 

- u 
(70) 



for < z < h, and 



u k (z) 



Rk 



(u - w ) rji coshrj 2 k(z + l 2 ) 
771 cosh r\ 2 kl 2 + r\ 2 co\hr\\kl\ svcfar\ 2 kl 2 



w 



(71) 



for — 1 2 < z < 0. The elastic (stray) energy per unit area 
of the film is now found by integrating over both parts 
of the film, since now a spontaneous strain exists in both 
of them: 

F s tray 16r)lT] 2 [l(uo - W ) 2 d 



A 



E 



[A1J1 + A 2 l 2 (wl/ul)] Tji+m 



1 - w /u 



14C(3) 



mm 



(74) 



in the same approximation as before, Trr) 2 l 2 /a •> 1, which 
enables us to replace coth by unity in (|7^). 

We find in the vicinity of T c2 , where r\ 2 = = 
^-2(A + 5A)/n < rji, A 2 < A x , w < u ), the equi- 
librium domain width 



'7r 3 AiZi 1 
14C(3)^ 



1 



x 



(75) 



We see that the period of the domain structure formally 
diverges when one is approaching T c2 from below but this 
behavior will be modified by the nonlinear effects. 

Let us check the behavior of the domain width at tem- 
peratures deep into the ferroelastic region for both parts 
of the film, where \A\ > SA. There rj 2 w 771, A 2 » Ai 
and 1 — wq/uq ~ <L4/2|A| <C 1, and we obtain 



2\A\ / 7T 3 

~SA\ 7V2C(3) 



d a t(h + h)- 



(76) 



We see that far below the temperature where a sponta- 
neous strain sets in the whole system, the period of the 
domain structure grows with respect to the Kittel period 
of initial domain structure qk, Eq. (|65|), as 



a 

a K 



2*/ 2 \A\ (h+l 



5A 



h 



1/2 



_ 9 3/2 



h + k 
h 



1/2 



T cl -T 



Tci - T c 



> 1. 



<-2 



(77) 



Since the period is linearly growing with lowering tem- 
perature, a oc and becomes very large, a 3> aK, and 
the replacement of the coth by unity becomes unjustified. 
The condition of applicability of Eq.(|77j) is the same as 
for Eq.(p5|) in the case of ferroelectrics. 



(2j + l) 3 Vi coth.r} 2 kjh + m coth-qikjh ' 



(72) 



IV. SUMMARY 



where 77* = \fM~Jji, i = 1,2, with M x = -2A, M 2 = 
-2(A + SA) and % = tt(2j + I) /a. Closer to T = T c2 
from below this expression becomes similar to that for 
the previous case, since wq — > 0. To find the total free 
energy one has to add the energy of the domain walls 



F dw _ yulAih + fiw^A 2 l 2 
A a 



(73) 



where A2 



A+8A 



1/2 



while Ai is given by (62) 



The equilibrium period of the structure is 



Summarizing, in an electroded ferroelectric or free fer- 
roelastic sample with a tiny inhomogeneity of either the 
critical temperature or temperature itself (i.e. in the 
presence of a slight temperature gradient and/or minute 
compositional inhomogeneity across the system) the do- 
main structure abruptly sets in when the spontaneous 
polarization appears in the softest part of the sample 
(i.e. the part with maximal T c ). This takes place when 
the difference in T c in the parts of the sample is just 
(1CP 3 — 10 _2 )K for the displacive systems, and even 
smaller, (l0~ 5 — 10 _4 )K, for the order-disorder systems. 
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The period of the structure then grows linearly with low- 
ering temperature and quickly becomes much larger than 
the corresponding Kittel period. 

This result does not depend on specific geometry as- 
sumed in the present model. Indeed, if local T c — T c (z) 
varies continuously it can be approximated by a piece- 
wise distribution of a sequence of "slices". Upon cool- 
ing the system first looses stability in the softest part of 
thickness l s , which is derived from the position of the 
boundary where local T c — 0, with respect to a domain 
structure with fine period oc \/T s . The domains extend 
into the bulk of the system and become wider with fur- 
ther cooling, since l s increases. In electroded sample 
there will be no domain branching and domain walls 
would run straight across all transformed slices. Oth- 
erwise, discontinuities would have resulted in very strong 
depolarizing field. If the overall inhomogeneity is small, 
the picture would obviously remain very similar to the 
two-slice model solved above. The same arguments re- 
main valid if the inhomogeneity were to have more com- 
plex form/distribution in a sample. The novel feature of 
the present effect in case of ferroelectics is that the depo- 
larizing field appears not due to surface charges, which 
are screened out by the electrodes, but because of the 
bulk inhomogeneity. In the case of ferroelastics, inho- 
mogeneity in the sample results in the bulk stresses that 
cause the splitting of the system into domains. In this 
case too the domain wall would run straight through the 
soft part of the crystal, since the discontinuities would 
result in large stray elastic stresses. 

We have shown that a very tiny temperature gradient, 
or a slight compositional inhomogeneity etc., would re- 
sult in practically any crystal eventually splitting into do- 
mains no matter how high the quality of it is. The rapid 
evolution of the domain pattern, found in the present pa- 
per, when it starts from very fine domains at T c , which 
then grow linearly with temperature to very lar ge s izes 
is similar to what have been reported in Ref. [|15| for 
~ 1mm thick TGS crystals. It would be very interesting 
to perform controlled experiments for the domain struc- 
ture close to the second order phase transitions. One 
could check, in particular, the basic assumption of the 
present theory that the electric fields (elastic stresses) ac- 
companying ferroelectric (ferroelastic) phase transitions 
even in slightly inhomogeneous media are compensated 
by formation of the domain structures rather than, for ex- 
ample, by screening of the electric field by charge carriers 
in ferroelectrics. Further understanding of the domain 
formation at phase transitions in real crystals is very im- 
portant given that many properties of ferroelectrics and 
ferroelastics, used in applications, are mainly determined 
by the domain structures. 
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